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ABSTRACT 

The evolution of Alfven turbulence due to three-wave interactions is discussed us- 
ing kinetic theory for a collisionless, thermal plasma. There are three low- frequency 
modes, analogous to the three modes of compressible MHD. When only Alfven waves 
are considered, the known anisotropy of turbulence in incompressible MHD theory 
is reproduced. Inclusion of a fast mode wave leads to separation of turbulence into 
two regimes: small wave numbers where three- wave processes involving a fast mode is 
dominant, and large wave numbers where the three Alfven wave process is dominant. 
Possible application of the anisotropic Alfven turbulence to the interstellar medium 
and dissipation of magnetic energy in magnetars is discussed. 

Key Virords: plasmas - turbulence - ISM: general - pulsars: general 



1 INTRODUCTION 



There is an extensive literature on the turbulence in the 
interstellar medium (ISM) that is based on incompressible 
MHD, in which case the only waves that can exist are Alfven 
waves. However, the ISM should be regarded as a compress- 
ible MHD medium at frequencies below an effective collision 
frequency, and as a collisionless, thermal, magnetized plasma 
at frequencies above the effective collision frequency. An im- 
portant qualitative change from the incompressible case is 
that there are then three modes, referred to as the Alfven, 
fast and slow modes in compressible MHD and as the Alfven, 
magnetosonic and ion sound modes in a collisionless ther- 
mal plasma. In principle, wave- wave interactions involving 
all three wave modes can contribute to the evolution of tur- 
bulence 

The conventional treatment of the evolution of tur- 
bulence in the ISM is in terms of a cascade of Alfven 
waves from a small wave number to a large wave num- 
ber through thr ee-wave int e ractio ns, as first proposed in - 
dependently by l lroshnikovl (^G^) and iKraichnanl il965l) . 
In the Iroshnikov-Kraichnan (IK) theory, the cascade pro- 
ceeds through collisions of oppositely directed Alfven wave 
packets, which leads to break-up of waves into smaller 
scales, corresponding to larger wave numbers. The IK the- 
ory predicts an energy spectrum E{k) ~ k~^^'^ in the in- 
compressible, isotropic MHD approximation. The isotropy 
assumption is valid only if the large scale magnetic field 
in the ISM is ignored. Inclusion of the large-scale mag- 
netic field changes the nature of wave-wave interactions 
and the spectrum of turbulence. An important effect is 



suppression of the cascade of Alfven waves along the 
mean magnetic field direction due to a kinematic restric- 
tion, required by the frequency and wave vector match- 
ing conditions llshebalin et alj Il983l: ISridhar fc GoldreichI 
1994h . A wave cascade perpendicular to the magnetic field 
is permitted, so that the turbulence evolves in a highly 
anisotropic manner, evolving much faster perpendicular to 
than parallel to the magnetic field, as sh own by numeri- 
cal MHD simulation JShebalin et alJll983D. This result was 
confirmed more recentiy~('Monteom^^ & Matthaeu^^M^ 
Ng fc Bhattachariec .1996,; G alti or ct al. .2000: Webb et alJ 
20001: iLithwick fc Goldreichll2003h . These studies lead to 
the conclusion that in the ideal, incompressible MHD ap- 
proximation, three-wave interactions of Alfven waves are al- 
lowed and lead to anisotropic turbulence with a spectrum 
E{k±) ~ kj^. However, there was some confusion in due to 
the three- wave matching conditi ons requiring that one mode 
must have a zero frequency ll^ i^a^^^oMrei^ ^94; 



Montgomery fc Matthaeu jll995l: iNg fc Bhattachariedll99' 



Galtier et al.ll2000l: IWebb et alJl2000l: ILithwick fc GoklreichI 
2003h . The assumption that a wave must have a nonzero fre- 
quency, combined with the three-wave matching condition 
requiring strictly zero frequency, led to the conclusion that 
Alfve n turbulence cannot evolve t hrough three- wave interac- 
tions JSridhar fc Go ldrcich 1994).^ However, this conclusion 
was later retracted ILithwick fc GoldreichI lioOSi) . with the 
zero-frequency, finite-fc^ disturbance being regarded as an 
acceptable component of the (perpendicular) turbulence. 

In this paper we consider the evolutions of Alfven turbu- 
lence in a collisionless, magnetized electron-proton plasma 
in the random phase approximation. As already noted. 
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the MHD approximation applies in the limit where the 
wave frequency is small compared with an effective colli- 
sion frequency, and it also requires that the wave num- 
ber be small compared with the inverse mean free path. 
These conditions are not satisfied, for example, for fluctua- 
tions on sufficiently small leng th scales in the ISM l)Ricketti 
ll990l:ISpang" er fc Gwinnlll99oli . Indeed, in essentially all as- 
trophysical and space plasmas of interest, waves at high 
wave number in the turbulent spectrum should be treated 
using kinetic theory. A more specific example of inter- 
est to us is the ultrarelativistic, superstrongly magnetized 
plasmas in magnetars, where magnetic energy dissipation 
through Alfven turbulence is thought to be important 
in transient X-ray emiss ion jThompson fc DuncanI Il995l: 
iThomspon fc Blaejll99^ . In such a medium, collisions play 
no role whatsoever, and the properties of the Alfven waves 
need to be treated using kinetic theory. 

An interesting formal problem arises when one ap- 
plies the theory of wave-waves interactions in a coUisionless 
plasma to the case where one of the modes has zero fre- 
quency. In the kinetic theory, the energy density in waves 
is proportional to their frequency, implying that there is 
no energy associated with the zero-frequency waves in the 
evolving perpendicular turbulence predicted by the MHD 
treatment. This conclusion cannot be correct, and the res- 
olution of the dilemma is not obvious. Possible resolution 
might be sought by identifying the perpendicular turbulence 
with some linear mode, but we are unable to identify any 
relevant mode. A more favorable approach is to modify the 
three-wave matching conditions by allowing for resonance 
broadening. This approach has already been developed in 
the context of ion sound turbulence where a similar problem 
arises: the three- wave matching condition cannot be satisfied 
for three ion sound waves, but such three wave interaction 
is possible when a broadening in frequency due to non linear 
effects is taken into account iTsvtovichlll97lL 119771) . The 
perpendicular turbulence is then regarded as a nonlinear 
quasi-mode of the medium, with the imaginary part of the 
frequency of the quasi-mode leading to the broadening of the 
resonance broadening. Any other form of nonlinear broad- 
ening that leads to a nonzero frequency would also overcome 
the problem in principle. Fortunately, in practice the zero- 
frequency wave does not lead to mathematical problems be- 
cause the evolution of the turbulence does not involve its 
frequency explicitly. 

Although we are concerned with the coUisionless regime, 
it is convenient to refer to the modes by their more familiar 
counterparts in compressible MHD. We denote the Alfven 
was as A, the magnetosonic (fast) mode as F, and the ion 
sound (slow) mode as S. We are aware of no treatments that 
include the S mode in three wave interactions involving the 
A and F modes, and we ignore this mode here. The simplest 
relevant approximation for the remaining A and F modes 
is the cold plasma limit, and most of our detailed results 
apply in this limit. The conventional cascade involves A —> 
A + A, and we concentrate initially on this interaction in 
the coUisionless case. However, even if the initial waves were 
purely Alfvenic (A), three- wave interactions, A —> A + F, 
A ^ F+F, would cause a mixture of modes to develop. As a 
result one also needs to consider the interactions F — > F+F, 
F ^ F + A, F A + A. The ratio of f to A in the resuhing 
turbulent spectrum is of interest, especially in view of the 



fact that the F can experience significant Landau damping, 
both by thermal particles and by cosmic rays. 

The properties of the Alfven mode and fast mode in 
a cold magnetized plasma are summarized in Sec. 2. Weak 
turbulence due to three-wave interactions involving three 
Alfven waves is discussed in Sec. 3. The effect of including 
the F mode is discussed in Sec. 4. Application to the ISM 
and to magnetars are discussed in Sec. 5 and 6, respectively. 



2 THE ALFVEN AND FAST MODES 

The properties of low-frequency waves in a coUisionless 
plasma are well known; the followin g summary of th ese prop- 
erties is based on the treatment bv lMelros3 (Il98d) . 

The general expression for the response tensor for a 
magnetized plasma with thermal distributions of electrons 
and ions, with thermal speeds Ve,Vi, involves modified 
Bessel functions with argument Xa = kj_Va/ila and plasma 
dispersion functions with argument Za = a;/\/2fc|| Va, where 
a = i,e denotes the species of particle with gyrofrequency 

= \qci\B /nia, and where where fcx,fc|| are the perpen- 
dicular and parallel wave vectors, respectively. The gen- 
eral form is approximated by assuming low frequencies, 
ui ^ Ljpi, and small gyroradii, ^ 1. The ion sound speed 
and the Alfven speed are identified as Vs — uipiVe/i^p and 
VA = {yii/LUpi)c. The nearly incompressible limit of MHD 
corresponds to the coUisionless regime Va va, which we 
assume to be satisfied. With these approximations the di- 
electric tensor reduces to the cold plasma approximation, 
which may be written in the form 

= ^J^^ [^'^ ~ + iv<.^e.jib^ , (2) 

where r]a = qa/\qa\ is the sign of the charge of species a, 
and b = B/_B is the magnetic field line direction, assumed 
to be the 3-axis. The low-frequency approximation corre- 
sponds to <C fii, where for simplicity we assume only one 
ionic species. The wave properties are then found by assum- 
ing the component \Kz2.\ is very large, and an expansion in 
uj/Q,i is performed. For present purposes, it is important to 
retain the corrections involving uj/Q,i only in the polariza- 
tion vectors. 

The dispersion relation, uj — (jj(k), the polarization vec- 
tor e, and the ratio R of the electric energy to the total wave 
energy for the Alfven mode (A) and the magnetoacoustic 
(fast) mode (F) may be approximated by 

uJA = \k\\\vo, GA = ei — cot^ 6*62, (3) 
lof = kvo, ~ 'q^. ^ei + iei (4) 
Ra = Rf = vl/2c, vl = v\/{l + v\/c), (5) 



ei = (cos sin 0), 62 = (— sin i/i, cos (/>, 0), (6) 



Anisotropic weak turbulence of Alfven waves in collisionless astrophysical plasmas 3 



with k = (fcx cos fcx sin (/), fc|| ) and tan 9 — k^/k\\. For 
v\ <^ c? one has ~ v\, and for sufficiently low den- 
sity or strong magnetic field one can have v\ ^ c? one 
has Vq ~ (? . Apart from the corrections proportional to 
uj'^ /Vll, these wave properties are then t he same as those 
be derived using relativis tic MHD theory (lAchterberslll983: 
iThomsDon fc Blaejll99l) . 

Thermal corrections to the dispersion relations and 
^ arise by retaining extra terms in the expansion in = 
k±Va/0,a and Za = Lu /\/2k\\Va- This is important when 
considering the damping of the waves. Landau damping by 
electrons dominates for the _F-mode waves, giving a damp- 
ing rate Tf oc Vs/vq- Landau damping of the F mode both 
by thermal electrons and by cosmic rays needs to be con- 
sidered. Landau damping for the A mode is much weaker, 
Ta oc (uj/fli)^, and can be ignored. 



3 WEAK ALFVEN TURBULENCE 

We apply the random phase formalism to three-wave inter- 
actions involving three A waves. 



3.1 Interaction of three Alfven waves 

A three-wave coupling between A waves may be denoted by 
A^ A' + A", where A, A', A" denote the three waves. The 
resonance or matching condition, which express conservation 
of energy and momentum from a semiclassical viewpoint, are 



k = k'-f-k", uj{k)=u;'{k')+u;"{k"), 



(7) 



with (jj(k) given by (|5Jl. These conditions then g ive /c|| — -j- 
fc'i' and I fey I = I fcfi I + 1 fcfi' I , which are satisfied only if one of the 
waves, say A" , has fc" = 0, corresponding to zero frequency. 
Then fcy = fcy implies that three-wave interactions do not 
affect the distribution in fcy and affect only the distribution 
in kx. 

The interpretation of the zero- frequency mode 
has been discussed extensively in the MHD con- 
text iMonteomerv fc Matthaeuslll995l: iBhattachariee fc Nd 



l200ll:lGaltier et alJl200ll2002l: iLithwick fc Goldreichll2003D 

A specific problem arises with the zero-frequency mode 
in the kinetic theory formalism because a zero-frequency 
disturbance carries no energy, seemingly precluding any 
transfer of energy in the turbulent cascade. The resolution 
of this dilemma that we propose here is that the nonlinear 
wave-wave interactions produce a natural spread in wave 
dispersion, i.e. growth or damping due to nonlinear inter- 
actions causes a broade ning in frequency in th e dispersion 
relations and @ jTsvtovichI Il97li IIQTtI) . The idea 
is that the frequency matching condition lO expresses 
energy conservation in the limit when the nonlinear 
transfer rate is negligible in comparison with the wave 
frequencies, and using it to infer a zero-frequency mode 
is inconsistent because wave frequencies smaller than the 
nonlinear transfer rate are not physically meaningful. An 
alternative version of this idea is as follows. For turbulence 
associated with waves in a specific mode, the turbulent 
energy is proportional to the Fourier transform of the 
correlation function of the wave amplitude, U(k,ijj) say. If 
the turbulence is sufficiently weak, such that the nonlinear 



effect on the dispersion relati on is neg ligible, then one has 
U{k,ij) = UCk)2nS{uj~uj(k)) jMelrosj tl986l. where (7(k) is 



the energy density per unit volume of k-space. The effect of 
the nonlinear transfer rate, F, is to broaden the J-function 
distribution into a Lorentzian profile centred at a frequency 
that differs from the usual dispersion relation, uiCk), by a 
nonli near frequency shift, with a width ~ |r| (iTsvtovichI 
'l97lV The inclusion of nonlinear broadening allows the 
three-wave matching condition to be satisfied by three 
A waves without requiring one of them to have a zero 
frequency. 

It is convenient to introduce the frequency mismatch 
between the three A waves: 



An = (|fc||| - |fc[|| - IfcJ'j) 



(8) 



For the third wave to have a nonzero frequency, due to a 
nonlinear correction or due to any other effect, Afl must 
be nonzero. A waves that all propagate either forward or 
backward cannot interact through the three-wave process. 
For example, if all three waves propagate forward fc|| > 0, 
fc'i > and fc'i' > 0, the mismatch would be zero, ASl = 0, 
but AQ 7^ is necessary when the nonlinear interaction is 
included, forbidding three forward-propagating waves from 
interacting. The inclusion of Afi 7^ 0, or uj" 7^ |fc||'|wo, is 
important conceptually, but it turns out to be unimpor- 
tant in treating the evolution of the A waves due to three- 
wave interaction. This is because the kinetic equation for the 
propagating Alfven waves does not depending explicitly on 
the frequency of this wave. One needs to appeal directly to 
the nonlinear broadening only when considering the kinetic 
equation for the perpendicular (zero-frequency) component 
of the turbulence. However, this is not considered in the 
MHD treatments, and we do not consider it here. 



3.2 Three- wave coupling 

In a semi-classical formalism (iTsvtovichI I1977I: iMehosj 
Il986h . the three- wave coupling is described by a probabil- 
ity that depends on the properties of the three waves and 
on the quadratic nonlinear response of the plasma. The 
quadratic response tensor, Qiji, is written down in the Ap- 
pendix. The probability can be written as Uaa'A" times a 
5-function that expresses the matching conditions Q, and 
is written below as (27r)''5''(fc — k' — k"). A detailed calcu- 
lation gives uaa'A" = '^^RaRa' RA"\ot^^, ^n?' /{el^^^'^"), 
with a^^/^/, = eAiCA' ^he projection of the non- 
linear response tensor onto the polarization vectors of the 
three waves. Some further details are discussed in the Ap- 
pendix. We find 



UAA'A" — 



k±c 

LU 



2 sin < 



COs{<l)" — (j)') + — - COS <j)' 



(9) 



where (p — is assumed, nip is the proton mass, assum- 
ing an electron-proton plasma, and rp = re{me/mp), with 
re = /'ineomeC^ ~ 2.6 x 10~^^ m the classical radius of 
the electron. As already noted, for a strictly zero-frequency 
mode Lj" = 0, @ implies a zero three-wave probability, but 
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there is a corresponding factor of l/ui" that effectively ap- 
pears when one writes down the kinetic equation for the A 
waves, allowing one to use Q without needing to specify 
the value of u" . 



trum in some limiting c ases using a method similar to that 
used bv lZakharovl il984l) . The kinetic equation llOII includes 
integrals over the azimuthal angles cj>' , which may be 
rewritten using 



hi -\- h i 



3.3 Kinetic equation for Alfven waves 

The evolution of the occupation number (classically, the 
wave action) A'^(k) of A-mode waves is determined by the 
kinetic equation 



dA^(k) 
dt 



dk' die" 



(27r)3 (27r)3 



UAA'A"{K'i^'X)NA"{-k") 



X [iV(k) - 7V(k')] (27r)''<5'*(fc - fc' - k") 



X [iV(k) - N{k'] 



dk", Ua" 



(27r)3 



o • J.' • J." ^0 

ismip smq) 



— COS0 



(10) 



On the left hand side of IIOII . the derivative is inter- 
preted according to d/dt = d/dt + -Vg ■ 9/9x + (duj/dx) ■ 
d/dk, where Vg = duj/dk is the group velocity. On the 
right hand side of (llOt the occupation number for the 
double-primed mode is NA"(k") — Ua" (k") /huj" , and for 
uj" —> this term dominates over other nonlinear terms 
o(^N(k)N(k') that are omitted here (Lithwick & Goldrei<3 
|20o3)- Note also the mentioned cancelation of uj" in the 
product MAA'A"(k,k',k")Ai'A"(k"). The equality in (ESJ is 
derived assuming Ri fcy and k^ ~ k± — k" . This approx- 
imation remains valid provided that the natural spread in 
frequency is much smaller than the wave frequency. Eg llOII 
is similar to an equation derived bv lGahier et alj 12003) (cf. 
their Eq [10]) using the incompressible MHD formalism. 

3.4 Natural broadening in frequency 

Nonlinear interactions lead both to a frequency shift and to 
wave growth or damping, giving a nonlinear correction to the 
dispersion relation, and introducing a natural broadening in 
frequency, respectively. The frequency shift and growth rate 
can be derived from 

Ra 



Sui — 



Re a 



Ra y N 

Ima , 



(11) 



(12) 



where all terms on the RHS are evaluated at the values from 
the dispersion relation, and with (k) = e* (k)e-i(k)ai^ de - 
termined by the quadratic response tensor jMelrosdlloi^) . 
As already noted, the evolution of C/^"(k") (the zero- 
frequency mode) requires F 7^ and the expression 11211 
generally leads to a nonzero growth or damping rate. 



2fcxfcj 



II 2 



On writing iV(k) = N{k±,k\^)/k± and UA"{k") 
Ua" {k'j_) /k'^, the kinetic equation reduces to 



(13) 
(14) 



diV(fcx,fc||) 
di 

X / dk', dk'. 



hi 

h 



k I k", 



1 - Mc 



1 " 

1 - 



N{k^,k^\)~-^N{k'^,k\\) 



1/2 



(15) 



where fi'^ and /i" are functions of k'j_/k±, k'\_/k± given by 
1131 1 a nd 114[l. and the t erms involving (fcn/fcx)"^ are ig- 
nored jGaltier et alJl2002^ . The integration range is limited 
by |cos0'] ^ 1 and |cos<^"| 1, implying k'^ > k± — k'J^, 
k'l < k± + k'^, and k'^ < k± + k'l. 

The turbulence spectrum is derived by a standard ar- 
gument based on a dimensional analysis. Assume that the 
source that drives the turbulence and the sink where the 
energy is dissipated are widely separated in k. A stationary 
solution with a power-law A'"(fcj .,fcii) ~ fc7" is s ought from 
list by a conformal transform iZakharovlll984r) . In the in- 
ertial range, where there is neither source nor sink, there is 
a constant flux of energy from smaller to larger k, and the 
spectrum is stationary, dN/dt — 0. It is convenient to intro- 
duce dimensionless variables ^1 — k'j^/k± and ^2 = k'j^/k±. 



,11 -0 



Assuming UA"{k'j_) 



, Eq HSU becomes 



X (1 - 



-l3-\, 



Cl 1 - 



1/2 



-l\^ia + f3-A 



(16) 



The second form of the integral is obtained from the first by 
writing = 1/Ci, ^2 = Mc(6,C2) ^ t^'c{€.'i,&) and 

1 - (6,6) ^ er'[l - m"c(6,^2)]. Half the sum of the 
two forms gives a third form 



dN 
It 



dCidC2 6~''~'ei ( 1 



X ( 1 - er''~M ( 1 - 6 



j 1 1 - 

a + /3-4 



/;2 



1/2 



(17) 



A stationary solution exists for a + l3 = 4. For a = f3, 
corresponding to a single power-law solution, the conditions 
for stationarity then gives a = P — 2, i.e. 



3.5 Turbulence spectrum 

Assuming that the three-wave interaction is the dominant 
process in the turbulence, one may derive the energy spec- 



iV(fcx,fc||) ~ kj 



(18) 



The s pectrum has the same form as for inc ompressible 
MHD iBhattachariee fc Ng|l200ll : I Galtier et al I bOQO, 2003 ; 
iLithwick fc Goldreichll200^ . 
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4 TURBULENCE INVOLVING THE F MODE 

The F (magnetosonic) mode may play an important role in 
the Alfven turbulence through the decay and coalescent pro- 
cesses F <-» A+A', A ^ F+A' . Such processes have no coun- 
terpart in incomp ressible MHD. Although, the se processes 
were discussed bv lLivishits fc TsvtovichiilQTOtl . the broad- 
ening effect was not considered. In this Section we include 
the broadening explicitly and discuss the relative importance 
of these processes to three- Alfven- wave interactions. 

Consider the case where a _F-mode wave interacts with 
two A waves, F ^ A + A' , that all propagate forward or all 
backward. It can be shown that when the frequency spread 
is included the three wave interaction is forbidden. Let {uj, k) 
and (oi', k'), {cu" , k") be the wave frequency and wave vector 
of the F-mode wave and the two A-mode waves, respectively. 
These satisfy the three- wave condition (|7J . Assuming AQ > 
to be a broadening to the three-wave frequency condition 
as the result of nonlinear interactions, for fey > and fey' > 0, 
one has 

fe = fefi -I- fe" -I- fc|| = fc|| -I- fcfi', (19) 



^"0 



which gives rise to 



(20) 



This equation is also valid when the two A waves both 
propagate backward, i.e. fey < and fey' < 0. Since the 
high frequency waves must satisfy uj ^ AQ ^ 0, the 
right-hand side is negative, forbidding F ^ A + A! with 
the two Alfven waves both propagating either forward or 
backw ard. An analogous resul t applies for compressible 
MHD jThomspon fc BlaejlToog) . 

It follows that F ^ A-\- A \& allowed only when the A 
waves are oppositely directed, fe||fe||' < 0. Assumin g fey > 
and fc'i' < 0, in a similar way to the derivation of 1201 . one 
finds 



fe 



A I h 



1 - 



An 



2lo" ) 



(21) 



where fe± ^ for Af2 < 2a;" is explicit. The general expres- 
sion for the three-wave interaction is rather cumbersome, 
and we consider only the special case where one of the two 
A waves has a low frequency, lj" <^ to, uj' . For <^ = 0, one 
has (cf. Appendix) 



2eoeLUpi / k±uju)" 



(22) 



The frequencies of all three waves are subject to nonlin- 
ear broadening (cf. Sect. 3.1), but such frequency broad- 
ening does not lead to a qualitative change in the nature 
of three-wave interactions. Although a low-frequency u" is 
assumed here, unlike interactions of three A waves consid- 
ered in Sect. 3, the three-wave condition does not requires 
w" — > 0. The probability can be approximated by 



+ 0"), (23) 



where 1211 1 is used to eliminate fex. In I23II . we assume 
uj" S> AO,, which is relevant here as Mfaa' oc uj""^ favors 
a moderately high uj" . In this approximation, the frequency 
spread plays no role in three-wave interactions. 



An important point is that three wave interactions can 
convert purely Alfvenic turbulence into a mixture of A and 
F modes. The process A F A' creates the _F-mode 
waves and the process F ^ A + A' removes them. In a 
steady state there is a mixture of F and A waves with the 
ratio determined by the ratio of the rates of the processes 
A ^ F + A' &nd A ^ A' + A". The process A ^ F + A' 
is described by the same probability as for F A -\- A'; 
one simply reverses the signs of the relevant ks and uses 
a;(— k) — — ci;(k). The probability 1231 continues to apply 
with a replacement </!>' ^ = 0. The ratio of _F to A waves 
in the turbulent spectrum can be estimated from the ra- 
tio of the interaction rate that generates (or destroys) the 
F-modes to the interaction rate for purely A-mode turbu- 
lence. These rates are ^afa' oc uafa'[NfO^') + A'"^/(k")] 
and 7aa'A" c< uaa'A"Na", respectively. Using (1231 and @, 
one finds 

Wf " 
Wa. 



Iafa' 



lAA'A" 



where Wf and Wa are the energy density of the F and 
A modes, respectively. One assumes that Wa' ~ Wa and 
va ^ c. We distinguish between u" '^'faa' in (HSJ from 
uj" uj'aa'a" ill ©• We define a critical frequency, Uc, that 
separates two regimes: for uj > ujc the rate of A ^ A' + A" 
is faster than the rate of A ^ F + A' , and for uj < Uc the 
rate of A — > F -|- A' is faster than that of the pure Alfvenic 
process A ^ A' -I- A". Using 1241 . one derives the critical 
frequency as 



1/3 



1 + 



Wf\ 

Wa) 



1/3 



(25) 



The inner and outer scales, corresponding to frequencies uji„ 
and ujout say, must satisfy ujin > oJc > i^out for this criti- 
cal frequency to be relevant. The turbulence in the regime 
UJ < LJc must involve F F' + F" as well as the above dis- 
cussed processes (involving both A and F). Unless a domi- 
nant process is singled out, in which case one may use the 
analytical approach similar to that discussed in Sec. 3.5, 
derivation of the spectrum in this regime needs a numerical 
approach and is not discussed here. 



5 THE ISM AND MAGNETARS 

In this Section we discuss possible application of the forego- 
ing results to turbulence in the ISM and to magnetic energy 
dissipation in magnetars. In both cases, we concentrate in 
the regime u <C O^. In the case of the ISM, this regime 
corresponds to a length scale well exceeding the ion inertial 
length scale VA/^i where wave damping through cyclotron 
resonance becomes important. 



5.1 Turbulence in a collisionless ISM 

The foregoing discussion leads to two qualitative features 
with implications for turbulence in the ISM. One feature is 
that the anisotropic features of Alfvenic turbulence found in 
treatments based on the incompressible MHD also apply to 
the collisionless case. The other feature is that when F mode 
waves are included, the ratio of the rates of the processes 
F ^ A + A' and A ^ _F + A' to the rate of the process 
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A yl' + A" is large at sufficiently low frequencies and 
small at sufficiently high frequencies. 

The presence of the F mode in the turbulent spectrum 
is relevant to scintillations of compact radio sources such a s 
pulsars and intraday variable sources (IDVs) llRickettll99Cll . 
Scintillations result from density fluctuations, which are as- 
sociated with the F mode but not with the A mode. Hence, 
if there is no F mode component there should be no scin- 
tillations. Our results suggest that the F mode should be 
present at frequencies lower than the critical frequency 



1/3 



where we assume Wf ~ Wa and u>" ~ Ari/2. The ISM 
has multi ple components with a ra nge of temperatures and 
densities IjMcKee fc Ostrikej ll977h. In McKee & Ostriker's 
three-component model, the ISM consists of a hot low den- 
sity component (the HIM), a cold neutral cloud component 
that is embedded in the former with a much smaller filling 
factor, and a warm coronal component (the WIM) that sur- 
rounds the neutral cloud. For the HIM, the typical density 
and temperature are ~ 3 x 10'^ m~'^ and T ~ 5 x K. 
The plasma frequency is LOpi = 73 {rii /3 x 10 ^ m"^)^^ s~^ . 
For a typical magnetic field B = 5 x 10"^° T dHeilejllQS'^ . 
the Alfven speed is va = {Q,i/LLipi)c ~ 2 x lO^ms"^ with 
f2i = 5 X 10~^ s~^. In such a plasma, the ratio of the thermal 
pressure (nefcsT) to the magnetic energy density {Wb ~ 
S^/2yLto) is generally less than unity, /3 — neksT/WB < 1. 
Assuming that the broadening is related to the ratio of 
the wave energy density to the magnetic energy density 
through AQ ~ uj{Wa/Wb) (e.g. Tsytovich 1977), one may 
estimate the critical frequency as Uc ~ 2u)pi{WA/WBy^^ ■ 
There is no reliable way to estimate the wave energy den- 
sity, though there must be an upper limit Wa/Wb ^ 1. 
For uic to be in the strong magnetic field regime {uic < fli), 



one must have Wa/Wb < {0.i/2u>p 



10 for the above- 



nominated parameters. Eq 1261 can be converted to a length 
scale Ic — vaI'^c for magnetic fluctuations in the par- 
allel direction (relative to the mean magnetic field). The 
lower limit to Ic is the ion inertial length li = va/^i'- 
I, ~ {vA/uJp^)(WA/WB)~'''''^ > ~ 4 X 10"^ m. So, the three 
A-wave process becomes dominant only in the region near 
the inner scale of the turbulence h < l/fc|| < Ic and the 
three- wave process involving F waves dominates in the short 
wavelength region l/fc|| ^ Ic- 

The above conclusion is valid for the WIM as well. 
The WIM has a relatively high density rn ~ 10^ m"^, 
corresponding to a plasma frequency ujpi = 1.3 x 

{rii/lO^ m~^y^^ , and a relatively low temperature 
T ~ 5 X 10^ K. One then estimates /c,wim/Zc,him ~ 
(ni,HiM/ni,wiM)^''^ = ^i,wiM/^i,HiM ~ 0.06, where the two 
components (HIM, WIM) of the ISM are labelled explicitly. 
Thus, the anisotropic Alfven turbulence may exist in both 
HIM and WIM, favorably in the scale of turbulence near h. 
It is emphasized here that this result is derived under the 
assumption of weak turbulence and that the possibility of 
strong turbulence in the ISM remains open. 



5.2 Alfven waves in magnetars 

Alfven turbulence is relevant to magnetic dissipation in 
magnetars, which are supercritically magnetized neutron 



stars (with B > B^r « 4.4 x 10^ T). Soft gamma-ray 
repeaters (SGRs) and AXPs are thought to be magnetars. 
Bursting X-ray emission is believed to be powered by mag- 
netic energy not rotational energy. The enormous magnetic 
stress beneath the surface causes cracks in the neutron 
star's solid crust producing elastic (shear) waves which 
couple to both the A mode and F mode that propagate in 
the magnetosphere iBlaes. Blandford. Goldreich. fc Madaul 
Il989il . These large amplitude waves trigger a cascade 
producing turbulence leading to plasma heating. The 
favorable region where the wave cascade occurs is the 
closed field line region where perturbations at the two 
opposite foot points of a closed field line launch A waves 
that meet head on. Since the shear waves underneath the 
surface propagate tangentially, the A waves excited by the 
shear waves may propagate approximately perpendicular 
to th e field lines (Blacs, Blan dford. Goldreich, & Mad<^ 
Il989tl . Assuming a sharp boundary, the shear waves 
can couple to both the A mode and F mode. Here 
we only consider the former. The typical frequency 
of A waves can be estimated fr om uiq ^ V/Ad = 
10^(K /10'^ ms~M(Ad/l m)~^ s~^ jThompson fc DuncanI 



I1995D . where Va — (fi/p)^^^ is the shear velocity, fi is 
the shear modulus, p is the crust density, and Ad is the 
length scale of the crack. The shear velocity is typically 
Vs ~ 10® ms"^ for fj, ~ 10^^ J m"^ and p ~ lO^^'kgm'^ 
In the strong magnetic field limit, one has va ^ c 
and the dispersion relation fcyc « uio- The latter gives 
fc|| ~ 3 X 10"^ m-^ 

For ideal MHD to be valid one needs the electron num- 
ber density in the plasma to be sufficiently high that it can 
support the required induction current that is dominant over 
the displacement current ~ uiSE / /iqc^ , where 5E ~ 5B/c 
is the electric field of the wave. Since the maximum cur- 
rent it can provide is risec, this condition can be written as 
i^p 3> {ujSD.ey^'^ , which corresponds to rie ^ 3 x 10^* m""^ 
for LU ~ 10® s"^ and SQ^ = eSB/rUe ~ 1.7 x 10^° s"^ (for 
SB ~ 10"^B with B = 10" T). As there can be a wide range 
of frequencies due to Ad not being well constrained, ideal 
MHD may not be applicable and the coUisionless regime may 
be the more appropriate. 

For the wave cascade to be relevant to plasma heating in 
magnetars, the decay time of A waves must be shorter than 
the flow time (~ the rotation period/27r) of the plasma in 
the magnetosphere. The propagation of A waves near the 
surface is oblique, with 



fex 
fell 



V 



3 X 10^ 



(27) 



for V ~ 10®ms~^ Jsiaes. Blandford. Goldreich. fc Madal 
119891 iThomspon fc Blaejll998h . The three-wave process in- 
volving a F wave is important at loq < uJc, where the critical 
can be estimated in a similar way to I26II . In the limit vo c, 
one finds 



fc|l 



5B_ 



Q.i ~ 10' 



(28) 



where Rq = 10^ m is the star's radius. Here we assume AO. ~ 
uj{5B/Bf, SB/B = 10'^ and fc±/fc|| = 300. The three A 
wave process starts to dominate at r ^ lO^iio for ujo — 
10® s~^ > Wc. In the following we consider the three A wave 
process only. 
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Since wave decay through three A wave interaction is 
preferentially in the perpendicular direction with the prob- 
ability cx the three- wave interaction would be sig- 
nificantly enhanced as fcx/fc|| increases. The decay time can 
be estimated from IIUII as 

where ALy is the longitudinal (along the magnetic field) size 
of the volume concerned. Since the decay process favours a 
large k± , the decay rate increases as the turbulence becomes 
anisotropic. The typical luminosity of transient X-rays (reg- 
ular bursts) from SGRs is about 10^* W. A giant burst may 
reach as high as lO^'^ — 10*" W. For example, the recent 
giant burst from SGR1806-20 on the 27t h December 2004 
relea sed about 10**" J in less than a second iBorko wski et all 
|20o3). Assume both high and low frequency waves are ex- 
cited with the energy density of the latter being Wa" ■ We 
also assume ALy ~ r, where r is the typical radial distance 
where wave cascade occurs, and that Wa" ~ 10^^ J/AV, 
where AV ~ r"^ is the volume. For k^\/k± ~ 300, and 
va/c ~ 7.2 X 10»(Bo/10"T)(n7l020m-3)-i/2(i?o/r)3/2, 
where Ro = 10* m is the star's radius. Bo is the magnetic 
field on the surface and rii is the proton number density, 
one finds ~ 1 s for a radius much smaller than the ra- 
dius of the light cylinder r ~ 3.9 x lO^iio < Rlc = 
cP/27r = 3 X 10*7?o(P/6s), where P is the magnet ar pe- 
riod. For Wa" ~ lO^^J/AV, which may be needed for a 
giant flare, the dissipation occurs much closer to the sur- 
face, r ~ 5 X 10^ -Ro. Wave cascade may occur much closer 
to the surface if one assumes a higher plasma density. 

One possible channel for dissipation is through Landau 
damping by thermal electrons. Since the frequency of the 
A" mode is nonzero, one expects the wave cascade would 
increase fcy and that the damping becomes effective when the 
parallel scale becomes sufficiently small l/fc|| ~ c/S7i ~ 1.8 m 
at a radial distance r ~ 3.9 x 10^ -Ro as compared to the 
initial l/fc|| ~ 300 m. Efficient damping would lead to rapid 
plasma heating and hence thermal X-ray emission. 



6 DISCUSSION AND CONCLUSIONS 

We study three-wave interactions of A waves in the random 
phase formalism in the collisionless limit. We consider weak 
turbulence so that energy and momentum of three inter- 
acting waves are conserved in each interaction. It is shown 
that similar to the recent result in the ideal, incompressible 
MHD, three- wave interactions of A waves can occur and lead 
to anisotropic turbulence in the collisionless approximation. 
In the incompressible formalism, three A waves can inter- 
act only when one of them has zero frequency as required 
by the three-wave resonance condition. Here we treat the 
zero-frequency mode as a low frequency limit lj" <ti uj ^ ui' 
but nonzero. The three-wave resonance condition can still 
be satisfied due to natural broadening as the result of non- 
linear interactions. Since each interaction results in only a 
small change in fcy , compared to a change in kj_, the turbu- 
lence becomes anisotropic in the sense that fiuctuations in 
the wave energy density are elongated in the perpendicular 
direction (in the k space), with spectrum ~ kj^. This result 
is similar to that derived in the MHD formalism. Because 



three-wave interactions also lead to a cascade in fc|| , dissipa- 
tion of wave energy is possible through Landau damping. 

Three-wave processes involving the F mode have a 
counterpart in compressible MHD but not in incompress- 
ible MHD. Rather than a single decay, A A' + A" , we 
also consider A — > F + A' and F ^ A + A' which create 
and destroy a -F-mode component in the turbulence. These 
processes involving the F mode are most important at the 
lowest frequencies, that is, near the outer scale of the turbu- 
lence. At higher frequencies three- wave interactions involv- 
ing three A waves become dominant, with the cross-over 
determined by I2(jll . Our resul t are consistent with the re- 
cent result bv lchandrarj ||200^ . who considered weak com- 
pressible MHD turbulence including the fast mode under 
the assumption of a constant density. We note that Cho & 
Lazarian's (2003) recent numerical simulation of turbulence 
in a compressible MHD shows that the F mode contributes 
negligibly to the energy transfer in the turbulence except 
near the outer scale region. However, their calculation is in 
the strong turbulence regime where the energy and momen- 
tum of three waves need not be conserved. 

We discuss two applications: Alfven turbulence in the 
ISM, and in the magnetospheres of magnetars. An im- 
portant implication for the ISM is the anisotropic nature 
of the turbulence which may be related to the observed 
anisotrop y in density flu ctuations, as suggeste d by several 
authors iGoldreich fc Sridha r 1995; Bhat tachariee fc Nd 
I2OOI). It is shown here that at a smaller wave number, 
three-wave processes involving a F wave is dominant and 
the process including A-mode waves only is important at a 
large wave number. Only the F mode involves density fluc- 
tuations that are required for scintillation. In contrast to the 
three A wave interaction, there is no preferential direction 
(in k) for wave cascades involving an -F-mode wave. The 
F mode can be dissipated through cascade to larger wave 
numbers and damped near the inner scale region. 

In the application to magnetars, the cascade time is 
found to be short so that the cascade can occur well in- 
side the radius of the light cylinder. Thus Alfven turbu- 
lence through three- wave interactions provides an important 
channel for dissipation. One possible way to dissipate the A 
waves is through Landau damping as they cascade to waves 
of short wavelength comparable to the inner scale of the 
turbulence. 

The third MHD-like mode, the ion sound or S mode, 
should be considered in any extension of the analysis in 
this paper. Although there have been some discussions 
of wea k compressible M HD turbulence involving the slow 
mode (lKuznetsovll200 1|) . we are aware of no detailed results 
that include the effect of the ion sound mode in magnetized 
turbulence in the collisionless limit. This mode is strongly 
(Landau) damped in a thermal plasma, and we speculate 
that its inclusion might lead to an effective damping of 
Alfven turbulence through the process A ^ A' + S, with 
the 5* mode rapidly damped. This might lead to an efficient 
dissipation of magnetic energy in the application to magne- 
tars, but we do not discuss this point further here. 
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APPENDIX A: THREE- WAVE PROBABILITY 

In the random phase approximation, three-wave interac- 
tions M ^ P + Q, where M, P and Q represent the 
three modes, satisfy the energy (frequency) and momen- 
tum (wave vector) conservation, which is also referred to 
as the three-wave resonance condition, i.e. k = k' -|- k" and 
Uj^,, (k) = LOp (k') +u)q (k") . Using a semi-classical formalism, 
a collection of waves in mode M can be described by their oc- 
cupation number NM{k) in the 6-dimensional phase space. 
Three-wave interactions are then described by the probabil- 
ity W MPoik, k' k") = UMPo(k, k\ k ")(2-R)^SHk - k' - k") 
with (|Tsvtovichlll977l : lMelroselll98(ih 

UAipQ{k,k ,k ) = — — ^ — — - 



egi^„(k)^^(k')^Q(k")| 
>^\o!MPQ{k,k',k")f 



(Al) 



where aj,,^pQ = e*^iepjeQsQiis, cvija is the quadratic re- 
sponse tensor, bm is the polarization of mode M, and 
k = {u, k). The 5-function contains the three- wave matching 
conditions. The quadratic response tensor in the cold plasma 
approximatio n has the foll owing symmetric (k' <-> k" and 
j ^ s) form jMelrosdll98^ 



E 



kr (a). Is {a), lis 



+— r^s'(t^ ' {uj ) + —rl ' {uj)tI/ {u- ^ 



7,// 7 // 



(C) 



(A2) 



(c) 



' zx 



0, 



with T^"' given by J^J. Since t._ 
one has a^fpQ{k,k' ,k") — for the Alfven or F modes 
propagating parallel or antiparallel to the mean magnetic 
field. In this special case, they do not interact with each 
other. 

It can be shown that in the expansion 
0{uj/Q,i,Lu' /Qi,uj" /Qi), only a^jl is important. Since 

Tzz ~ 0{1), a^^^l is nonzero only for i = j = s — z. 
However, since caz ~ 0(1/^1^) x 0{v'i/v\) (caz = for a 
cold plasma), one has bazSa' z&a" z'^^^'^^^ ~ 0{l/Q,%). The 
nonzero components of a^^l correspond to one or two of 
the three indices being z. Then they produce terms of 
0(1/Ql) X 0{v'^Jv\) and 0{1/Ql) in aAA'A"- One can 
show that a'-^jl = for </> = 0. To the order 0{1/D.l), the 
relevant components are a^xx, cixxy, cixyx, O-yxx, Q-xyy, 
ctyxy, and ctyyx- Since rri^Ql 3> mpf^f, one may ignore the 
electron terms: 



E 



qgng 1 



E 



qiUi 1 



+ 0{mjmi). 



(A3) 



In the cold plasma approximation, the polarization vec- 
tors are effectively perpendicular to B and then all the com- 
ponents of aiji with any of the three indices equal to z gives 
zero contribution to a^j^-pg (fc, fc', fc"). The remaining terms 



give = -{e^n,/mlcn^)&];> with 



(n) 
ijs 



OLxxx ~ i {kyLJ + 2k'yu") CLO, 



^xxy 


:i i (cj" — (x)) k"cuj, 


~ (3) 
^xyx 


■ i 1 1 ll\ 

n I y — k,j.cuLU ) , 


yxx 


i i [ — kxCLuuj") , 


^xyy 


i ikyCLuu" , 


"■yxy 




^yyx 


^ ''' 1 "\ 1 " 


^yyy 


- {-k'^uj - 2k'^uj") cuj, 


-(4) ^ 

^xxx 


^0, 


^xxy 




-(4) 
^xyx 


s 2 [— ojfcy + Lu" {ky — fcy)] CLO 


~(4) 
xyy 


- 2 [uk'l + u;"(k'^ - fc")] cuj'^ . 


We only retain terms up to 0{uj" /lo 
in the expansion give 



(A5) 
(A6) 
(A7) 
(A8) 

(A9) 
(AlO) 

(All) 

(A12) 
(A13) 
(A14) 
(A15) 



- e'^'e^"' e(°>*a<*' -he(''e<°' e'°'*a(^' 

^ ^Ai'^A'j'^A"s^ijs ^ ^Ai^A'j^A"a'-^ijs 

^'^Ai^A'i^A"s'-^ijs ^ ^Ai'^A'j^A"s"-ija 
^ ^A'j^A"s'^xis '^"^ ^^A'j^A'is'^yjs 



^'^A'j'^A"s"a;j3 A' A" s'-^xj s ^ 



(A16) 



where the polarization vector is written into the form e ~ 
e'"' + e*-^' The second equality is derived for = 0. For 
fc|| ~ fc'i one has ] cose"] = |fc||'|/\/fcf|'^ + k'{^ < 1. Thus we 
have 



111 



kf + k'!' 



0{lo"^). (A17) 



The first term in Eq. 1A16II is dominant, which is used to 
calculate @. Here one assumes a single ionic component 
i — p (protons). 

The quadratic response for F ^ A + A' can be written 
as eFieAjeA'sCtijs ~ epieAjeA'sOifj's- One may derive 
using (A4)-(A11). 
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